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Abstract. The goal of this paper is to estimate directly the rotation 
and translation between two stereoscopic images with the help of five 
homologous points. The methodology presented does not mix the ro- 
tation and translation parameters, which is comparably an important 
advantage over the methods using the well-known essential matrix. This 

(f) ■ results in correct behavior and accuracy for situations otherwise known 

^ ' as quite unfavorable, such as planar scenes, or panoramic sets of im- 

C**** , ages (with a null base length), while providing quite comparable results 

\l ■ for more "standard" cases. The resolution of the algebraic polynomials 

resulting from the modeling of the coplanarity constraint is made with 
the help of powerful algebraic solver tools (the Grobner bases and the 

£ — ■ Rational Univariate Representation). 

o 
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O ■ 1 Introduction 
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The determination of the relative orientation between two cameras with the help 
of homologous points is the basis of many applications in the domains of pho- 
jrt ■ togrammetry and computer vision. The configuration often called "minimal case 

problem" takes the intrinsic parameters (i. e. the elements of calibration) of the 
camera as a priori known. Then only five points homologous are necessary to 
estimate the remaining three unknowns of rotation and two ones of translation 
(up to a scale factor). This problem has been dealt by many authors, and most 
of recent methods published provide a resolution based on the properties of the 
essential matrix. Even if its use simplifies remarkably the problem of the relative 
orientation, in some cases, due to the fact that all parameters of rotation and 
translation are mixed, this is the origin of geometric ambiguousnesses. So as to 
improve this point, we propose in this article a model that separates completely 
the rotation and the translation unknowns. We show that the major advan- 
tage of this model is that it allows to solve degenerate problems such as pure 
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rotations (null translation). We use an algebraic modeling for the coplanarity 
constraint, through a system of polynomial equations. We solve them with the 
help of powerful algebraic solver tools, the Grobner bases and the Rational Uni- 
variate Representation. So as to assess this new approach, three cases have been 
processed: a classical case, a planar scene, and a case where the base length is 
close to zero. We will see that the new method is still accurate even for the last 
two cases - quite unfavorable - configurations. We will also compare with the 
Stewenius's algorithm and see that in planar scenes the new algorithm is more 
accurate. An evaluation on real scenes will finally be presented. 

2 Historical background of the five points relative pose 
problem. 

It was for the first time demonstrated by Kruppa [1] in 1913 that the direct 
resolution of the relative orientation from 5 points in general contained at most 11 
solutions. The described method consisted to find all intersections of two curves 
of degree 6. Unfortunately, one century ago, this method could not lead to a 
numerical implementation. Lately in [2], [3], [4], [5] it has been demonstrated that 
the number of solutions is in general equal to 10, including the complex solutions. 
Triggs [6] has provided a detailed version for a numeric implementation. Philip 
[7] presented in 1996 a solution using a polynomial of degree 13, and has proposed 
a numeric method to solve his system. The roots of his polynomial give directly 
the relative orientation. Philip has exploited the constraints on essential matrix. 
Philip's ideas have been followed in 2004 by Nister [8] who has refined this 
algorithm, has obtained a 10*' 1 order polynomial and has given a numerical 
resolution using a Gauss- Jordan elimination. Since then, number of papers tried 
to give some improvements to the method of Nister, notably Stewenius [9] that 
has provided a polynomial resolution using the Grobner bases. Many papers 
have proposed some modifications to the method of Nister in view of a numeric 
improvement [12], [13], or for a simplification of implementation [10], [11]. 

3 Geometry review of Relative Orientation 

In this section we recall the various ways to present the geometry of relative 
orientation, that consist in the determination of the translation and the relative 
rotation between two images of a scene having a common informational part. 
In general, the position of the first camera is taken as the origin of the system 
Si (Fig. 1) and therefore the position of the second camera (S2) is calculated in 
relationship to the first one. O named as the principal point of the camera, / is 
the focal length. A is the world point, and the image projection of A on the left 
image is a with coordinates (x a ,y a , f) T , and a' (x a ',y a i , f) T in the right image. 
The vector of translation T (T x ,T y ,T z ) is the basis that relies the optic centers 
of the cameras (Si and S2). R is the relative rotation between the two cameras. 
A way for modelling the relative orientation is known as condition of coplanarity. 
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Fig. 1. Geometry of relative orientation 



This constraint has been often used by the community of computer vision since 
three decades. As pictured in the Fig. 1, the condition of coplanarity between 
two images expresses the fact that the vector V\ , the vector V 2 (expressed in the 
reference of V\), and the vector of the translation T are in the same plane, called 
the epipolar plane. One can translate this condition by a null value for the triple 
product between these 3 vectors. In other words: 



V 2 -{RVi/\ T) = 



(1) 



4 Algebraic modelling of the five-points problem 

In this section, different ways for algebraic modelling of the relative orientation 
are recalled. The goal is to obtain a polynomial system, so as to use the powerful 
mathematical tools developed for solving such systems. 

The coplanarity constraint (equation 1) under its algebraic shape is expressed 
by the equation: 
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In this equation the unknowns are the matrix of rotation R and the translation 
T. Different ways exist to parameter the system so as to obtain polynomials 
with the rotation and the translation as unknowns. In the present part the main 
modelling solutions for the rotation and the translation to be used in this research 
are described. 

4.1 Modelling of the translation 

The translation of unity length between the two centres of the cameras may be 
understood as imaging on the unity sphere its center. The translation has only 
2 degree of freedom, and for that reason, with the relative orientation, the scale 
cannot be determined. The equation of the unity sphere is the following : 



Tl + Tl 



■ Ti = 1. 



(3) 
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The advantage of this constraint of normality on the translation is that it is 
quite possible to work even with a very small translation, allowing to compute 
precisely the rotation when the translation is null. In this case, as the base T 
is null, the two homologous vectors V\ and Vi are deduced from each other by 
a rotation R, so that RV X A^ = "0. Thus the triple product T • (RVi A V 2 ) 
is null whatever T . The fact that we force the translation to be unity prevents 
that T = and therefore to suffer numeric instabilities. This implies in turn 
that the rotation will be correctly estimated in any case. 



4.2 Modelling of the Rotation in 3D space 

The rotation matrix (R) in the 3D space has 3 degree of freedom. It is thus pos- 
sible to express it with 3 parameters. However several representations with more 
than 3 parameters exist. The algebraic model will depend on the chosen repre- 
sentation. In the following part the main models for the coplanarity constraint 
are described. 



Representation of the rotation using quaternions (4 parameters) A 

quaternion is composed of 4 parameters, q = (a, b, c, d)', with the vector part 
being (b,c,d). The quaternions provide a simple representation of the rotation. 
Indeed with the parameters of a unity quaternion , the matrix of rotation can 
be expressed in the following manner : 



1 - 2 c 2 - 2 d 2 2 be - 2 da 2bd + 2ca 
2 be + 2 da 1 -2b 2 -2d 2 2cd-2ba 
2bd-2ca 2cd + 2ba 1 - 2 b 2 - 2 c 2 



(4) 



With such a model the number of unknowns for the rotation is 4. 



Representation using Thompson rotation. Another efficient way to repre- 
sent the rotation with three parameters is given in Thompson's paper [14]. 
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where A = 1 + ±(A 2 + A t 2 + v 2 ) and A' = 1 - \{X 2 ■ 

unknowns for the rotation also res 
rotation matrix exist, such as the Cayley transfom, often used in robotics [15]. 



r fi 2 + iy' 2 ) and A' = 1 - |(A 2 +/i 2 + v 2 ). With such a model 
the number of unknowns for the rotation also resumes to three. Other models of 
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4.3 Algebraic modelling of the coplanarity constraint 

Model with 7 equations With the model of rotation matrix given by equation 
4, every couple of homologous points (a and a') gives an equation of the type : 

(xa' l (-T z (2bc+2da) + T y (2bd-2ca)) + ya' l (T z (l-2c 2 -2d 2 )-T x (2bd-2ca))+ 
za t (-T y (l-2c 2 -2d 2 )+T x (2bc+2da)))xa2+{xa l (-T z (l-2b 2 -2d 2 )+T y (2cd+2ba))- 
ya',(T z (2bc-2da)-T x (2cd+2ba)) + za' t (-T y {2bc-2da)+T x (l-2b 2 -2d 2 )))ija 2 + 
{xa! l (-T z (2cd-2ba) + T y (l-2b 2 -2c 2 ))+ya' l (T z (2bd+2ca) -T x (l-2b 2 -2c 2 ))+ 

za' t (-T y (2bd + 2ca) + T x (2cd - 2ba)))za 2 = (6) 

5 couples of points gives 5 equations. It is otherwise necessary to add along with 
this equations the constraint of normality on the quaternion, as well as on the 
translation. 

T 2 +T 2 +T 2 = l, (7) 

a 2 + b 2 +c 2 +d 2 = l. (8) 

Finally, the unknowns of the system are [a, b, c, d, T x ,T y , T z ] 

Model with 6 equations While using the Thompson rotation matrix, the 
rotation is expressed with 3 parameters. The system will have 6 unknowns, 
considering the three parameters of translation. The polynomial expressing the 
coplanarity constraint for a couple of homologous points, taking for model the 
Thompson rotation, is the following : 

(xa t (-4T z v - 2T 2 A/i - 47> + 2T y Xv) + ya,(AT z + T z \ 2 - T zf i 2 - T z v 2 + 

4T x /z - 2T x \v) + zai(-4T v - T y \ 2 + T y fi 2 + T y v 2 + AT x v + 2T x X^))xa' \+ 

(xai(-4T z + T z \ 2 - T z ^i 2 +T z v 2 + AT y \ + 2T y [iv) + ya t (-4T z i> + 2T z Apt - AT X X- 

2T x fj,v) + za,i(AT y v -2T y Xfi + 4T x -T x \ 2 +T x fi 2 -T x v 2 ))ya' , + (xa l (4T z X~2T z fiiy 

+ AT y - T y X 2 - T y[ ? + T y v 2 ) + j/^ (4T^ + 2T z Xv - AT X + T X X 2 + T x[ i 2 - T x v 2 )+ 

za t (-4:T y fi - 2T y Xv - 4T X X + 2T x fiis))za' t = (9) 

The constraint of normality on the translation (equation 3) is added to these 5 
equations. So the system has 6 equations and 6 unknowns [A, /i, v, T x , T y ,T z }. 

In conclusion of this section, we have built two polynomial systems, where 
the translation and rotation parameters are distinct and correspond to separated 
unknowns. Next, we show how to solve this type of polynomial systems. 

5 Resolution of the polynomial systems 

The ways to solve the polynomial systems are widely published [16], [17], and 
are briefly recalled for the reader not familiar with this topic. The resolution 
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of a polynomial system consists in finding the zeros of an algebraic equation 
system such as : P{x) = with P = [p\,p 2 , ■■,Pn) where pi is a I — variable 
polynomial x — (x±,x 2 , ■ ■-,xi) over the field C of complex numbers. Diffcrents 
types of solvers for polynomial equations exist, such as analytic solvers, subdi- 
vision solvers, geometric solvers, homotopic solvers and algebraic solvers [18]. In 
this paper the focus is on algebraic solvers, that exploit the known relationships 
between the unknowns. They subdivide the problem of the resolution into two 
steps : the first consists in transforming the system into one or several equivalent 
systems, but better adapted, and this constitutes what one will call an algebraic 
solution. The second step consists, in the case where one works in one subfield 
of the complex field, to calculate the numeric values of the solutions from the 
algebraic solution. We will sec now briefly the principal tools used in this paper 
for solving polynomial systems. But first, some remainders of geometric algebra 
are necessary. 

5.1 Notations 

Q [Xi, X 2 , ..., X n ] is the polynomial rings with rational coefficients and unknowns 
X u X 2 ,...,X n . S = P U P 2 ,...P 3 is any subset of Q [X U X 2 , ...,X n }. A point 
x G C" is a zero of S if Pi(x) = Vi = 1, 2, ..., s. The variety of P is the set of 
all common complex zeros : 

V(P) = {(oi, ..., o„) G C" : pi(ai, ..., a n ) =0foralll<i< s}. (10) 

The ideal 1 generated by a finite set of multivariate polynomials < P 1 ,P 2 , ...,P S > 
is defined as: 

n 

l={J2h i P i \h i €Q[X 1 ,X 2 ,...,X n }}. (11) 

The ideal contains all polynomials which can be generated as an algebraic 
combination of its generators. An ideal can be generated by many different sets 
of generators, which all have the same solutions. 



5.2 Construction of the algebraic solver : an introduction to the 
Grobner bases. 

A Grobner basis is a set of multivariate polynomials that has " nice" algorithmic 
properties. Every set of polynomials can be transformed into a Grobner basis. 
This process generalizes three familiar techniques : the Gauss elimination for 
solving linear systems of equations, the Euclidean algorithm for computing the 
greatest common divisor of two univariate polynomials, and the Simplex Algo- 
rithm for linear programming [3] . The Grobner bases were developed initially by 
B. Buchbcrgcr in the years 1960 [19]. The first step, when we want to compute a 
Grobner basis, is to define an monomial order. For polynomial rings with sever- 
able variables, there are many possible choices of monomial orders. Some of the 
most important orders are given in the following definitions. 
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Definition 1 Lexicographic Order (Lex) 

X? 1 ... X%» < Lex Xf 1 ... Xt <* 3 t Q < n 



&i = Pi, for i = 1,2, ...,i - 



Pk 



Definition 2 Degree Reverse Lexicographic Order (DRL) 

xp...xz» < DRL x*...xfr & xp k x^...x-^ < Lex xpx-' 3 

The DRL ordering, although somewhat nonintuitive, has some interesting com- 
putational properties. 

The following terms and notation are present in the literature of Grobner bases 
and will be useful later on. The degree of a polynomial P, denoted DEG(f), is 
the highest degree of the terms in P. The leading term of P. denoted LT(P), is 
the term with the highest degree. The leading coefficient of P denoted LC(P) is 
the coefficient of the leading term in P. Finally Grobner bases can be defined: 

Definition 3 Fix a monomial order > on Q [Xi, X2, ..., X n ], and 
let X C Q [Xi, X2, ..., X n ] be an ideal. A Grobner base for X (with respect to >) 
is a finite collection of polynomials G = {51, ...,gt} C X with the property that 
for every nonzero /el, LT(f) can be divided by LT(gi) for some i. 

Two principal questions immediately arise from this definition : 

1 . the existence of Grobner bases for any ideal X. 
Hilbcrt's Basis Theorem says that : Every ideal X has a Grobner basis G. 
Furthermore, X =< <ji, ..., g t >. 

2. the issue of uniqueness of the Grobner bases. 
The Buchbcrgcr theorem proves that if we fix a term order, then every non- 
zero ideal X has an unique reduced Grobner basis with respect to this term 
order. 

There are several possible algorithms to effectively compute Grobner bases. The 
traditional one is Buchbergers algorithm, it has several variants and it is im- 
plemented in most general computer algebra systems like Maple, Mathematica, 
Singular [20], Macaulay2 [21], C0C0A [22] and the Salsa Software [23]. In this 
paper we use the Salsa Software with the F4 algorithm [24]. The Faugere F4 
algorithm is based on the intensive use of linear algebra methods. 

5.3 Application of Grobner bases for systems solving 

Grobner bases (G) give important informations about the initial system of poly- 
nomial equations. Here some of these informations with their application on the 
2 systems of polynomials defined in Section 4.2 and Section 4.3 are summarized: 

1. Solvability of the polynomial system. It is easy to check if it is possible to 
solve the system of polynomials or not with the help of the Grobner bases G. 
If G = {1}, the system has no solution. We check this on our two systems, 
and we find that: G ^ {1}. In other terms V is not empty. 
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2. Finite solvability of polynomial equations. It is easy to see whether the system 
has a finite number of complex solutions or not : we just check that for each 
i, 1 < i < n, there is an m; > such that x™" 1 = LT(g) for some g £ G. 
This type of system is called zero-dimensional system. In this case, the set 
of solutions does not depend on the chosen algebraically closed field. If we 
apply this on two systems, we find that the dimension of the two systems is 
zero. So the set of solutions is finite. 

3. Counting number of finite solutions of the polynomial system. One important 
information is that the Grobner basis also gives the number of solutions of 
the system. If we suppose that the system of polynomial equations P has a 
finite number of solutions, then the number of solutions is equal to the car- 
dinality of the set of monomials that are not multiple of the leading terms 
of the polynomials in the Grobner basis (any term ordering may be chosen) . 
This monomials are called basis monomials or standard monomials (B). 
For the system of 7 polynomial equations (Section 4.2), the number of the 
standard bases is equal to 80. It means that the system has a total of 
80 complex solutions. This high figure is due to the symmetries of the 2 
equations of constraints, for example the quaternion Q(a,b, c, d) is equal to 
Q'(—a, —6, — c, — d). It has been reminded previously in the literature that 
the maximum number of solutions for the essential matrix is 10. Every es- 
sential matrix gives 4 families of solutions, rotation and translation, which 
provides 40 families of solutions. In the present system, with the symmetry 
on the quaternion, there are twice more solutions, therefore the total is 80. 
Using the system of polynomial equations defined in Section 4.3, the stan- 
dards bases of this system are the following (in the DRL order): 

B = [1, T Zl T y ,T x , v, /i, A, T 2 , T y T Zl Ty, T x T z ,T x T y ,uT z ,pT y , vT x , v 2 \ {iT z , 

IxTy, fj,T x , uv, /i 2 , \T Z , XT y , XT X , \v, A/i, A 2 , T%, T y T z , T 2 T Z , T X T Z , 

T x T y T z ,vT z , vT y T z , is 2 T z ,fiT z , fiT y T z , fivT z , AT 2 , \vT z ] 

(12) 

Which makes a total of 40 bases and therefore 40 solutions. In the present pa- 
per the Salsa library has been used. It is quite possible to use the resolutions 
according to Stewenius [9] and [25], but it requires a heavier implementation 
work. 

5.4 Finding the real roots of the polynomial systems 

Once the Grobner basis is calculated, different ways exist to find the roots of the 
system of polynomial equations, e.g. the method that solves the polynomial sys- 
tems with the help of elimination and lex Grobner basis. Another most popular 
way to solve polynomial systems is via eigenvalues and eigenvectors, with the 
help of standard monomials [16]. Here the emphasis is put on the other method, 
called the Rational Univariate Representation (abbreviated by RUR) [26]. Rep- 
resenting the roots of a system of polynomial equations in the RUR was first 
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introduced by Leopold Kronccker [27], but started to be used in computer alge- 
bra only recently [28], [29], [26]. The RUR is the simplest way for representing 
symbolically the roots of a zero-dimensional system without loosing information 
(multiplicities or real roots) since one can get all the information on the roots of 
the system by solving univariate polynomials. Let P(X) = (< P\, P%, ...P s >) 
where Pi £ Q [X\, X 2l ...,X n ] be a zero-dimensional system with its solution 
set V = P _1 (0), Rational Univariate Representation of V consists in expressing 
all the coordinates as functions of the roots of a univariate polynomial such as : 

fo(T) = 0,X 1 = ll^l,X 2 = ^l,....,X n =^l (13) 

where fo, fi, f 2 , ■■■, fn,q S Q [T] (T is a new variable). Computing a RUR reduces 
the resolution of a zero-dimensional system to solving one polynomial with one 
variable (f t ) and to evaluate n rational fractions ( ( T -! , i = 1, •■■, n) as its roots. 
The goal is to compute all the real roots of the system (and only the real roots), 
providing a numerical approximation with an arbitrary precision (set by the 
user) of the coordinates. Many efficients algorithms have been implemented to 
calcultate RUR. More details are easily found in the literature, but a complete 
explication can be found in [30], [26], [31]. An implementation of the Rouillicr 
algorithm for RUR computation can be found in the SALSA software [23]. 



6 Algorithm Outlines 

Now, the different steps of our algorithm for the calculation of the relative ori- 
entation are described. 

Step 1 : 5 couples of homologous points are randomly selected with the RANSAC 
method [32] [33] [34]. 

Step 2: Build the system of polynomial equations. 

Step 3: Solve polynomials system. In the present paper the Salsa library has 
been used (it is quite possible to use the resolutions according to Stewenius [9] 
and [25] , it requires a heavier implementation work, but has the advantage that 
a great part of the calculation is made off-line) . 

Step 4 '■ Identify the solution with a physical sense. The ambiguity resolution 
may be done through the use of a third image [8], but we prefer to be able to 
work with only two images. It is important to find the "true" solution in this 
very large set, and it is necessary to inject information bound to the geometry 
of the scene. We proceed in this way : 

— when intersecting the rays relative to all homologous couples of points, one 
keeps the solution where the rays cut themselves in front of the image, 

— for the 5 randomly selected couples, we calculate the distance to the world 
points. We keep only the solutions that give a depth superior to the value of 
the normalised baseline, i. e. 1, the other ones being considered as unrealistic, 

— the last step consists in selecting the solution which fits with the highest 
number of points. This hypothesis requires of course to have more than 
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five points. Other methods exist to find the good solution among all those 
produced by the direct resolutions, but in general they consist in using a 
third image [8]. 

As an example, Fig. 2 illustrates the different steps. The reference image is in 
blue frame, all real solutions are presented in green frame, and the final solution 
has a red frame. 




Fig. 2. Representation of all real solutions 



7 Results and Evaluation 

Here we present the results of an experimentation on both synthetic and real 
data. 



7.1 Experimentation on Synthetic Data 

To quantify the performances of the presented method, synthetic data have been 
simulated. The parameters used for the simulations, arc the same as Nister's 
ones. The images size is 352 x 288 pixels (CIF). The field of view is 45 degrees 
wide. The distance to the scene is equal to 1. Several cases have been treated : 

1. Simple configuration : the baseline between the 2 images has a length of 0.3, 
the depth varies from to 2. 
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2. Planar Structure and short baseline (0.1) : a degenerate case where all sim- 
ulated points are on the plane Z = 2. 

3. Zero translation : the configuration of the points is the same as in the simple 
configuration, the main difference is that the baseline length is null. 

In each configuration a Gaussian noise with a standard deviation varying between 
and 1 pixel is added. The results are average of 100 times indcpcndant experi- 
ments. For each situation the minimal case only has been treated, corresponding 
to the minimum number of points required (5). No least square adjustment has 
been done. The geometry of the different configurations is illustrated in the Fig. 
3. 




(a) (b) (c) 

Fig. 3. (a) Easy condition (b) Planar Condition (c) Zero translational condition 



Results in the so-called easy configuration . As a comparison our method 
has been confronted with Stewenius's one, thanks to the codes that he kindly 
made downloadable on his website [35] . This allowed us to use it as a reference for 
the algorithm of 5 points. Two sorts of translations have been treated, one in X 
(sideway motion) and one in Z (forward motion). Our results are mostly similar, 
even slightly better. Remark: In these simulations, it is important to specify that 
the rotation between the two optical axes is always very well determined. The 
difference between the different methods is the precision of evaluation on the 
orientation of the base, so this is the assessment that we have used. 



Planar Structure and short base . Several surfaces are known as "dan- 
gerous" [36] the reason of this appellation is due to the fact that if the points 
chosen for the evaluation of the relative orientation are on this kind of surface, 
the configuration becomes degenerate. In the following, one of the most unfavor- 
able configurations has been chosen. We note that the method of the 5 points 
of Stewenius is not robust in the sideways motion case. Besides, Sarkis [13] has 
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(a) 



(b) 



Fig. 4. Error on the base orientation (in degree). Easy Case, a) sideway motion, b) 
forward motion. 



shown this weakness of the algorithm, and concluded that it is better in such 
cases to use an homography. On the other hand, with the method presented here 
this kind of configuration does not lead to any problem. 



§ °° 

1 ■» 

I 30 

B 20 
10 


<-v 


^-"""^ 




\ 


I -a- New 5pt 
-*- STEW06 









1 3 

S. 5 




(a) 



(b) 



(c) 



Fig. 5. a) Error on the base orientation (in degree) and standard deviation, planar case, 
sideway motion b) idem, for a forward motion, c) Error on the relative orientation, non 
planar case, null base length. 



Results for a null base length Even with a null translation, the rotation 
is very well defined. This is probably due to the fact that the parameters to 
estimate in our initial equations are completely separated. The rotation is not 
mixed with the translation. This may explain why a translation of zero length 
does not affect the results. Fig. 5(c). 



7.2 Tests on Real Images 

So as to test the algorithm presented in section 6, we have used the recently 
available image base set up by ISPRS [37] for relative orientation tests. We 
have checked especially our ambiguity resolution so as to find the good physical 
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solution. The mean error on the baseline orientation is equal to 5.25° and mean 
error on the rotation is 1.26°. At the present time, an average of 20 seconds on a 
1.65 GHz machine is necessary to get the solution. This time of calculation may 
be to improved by using offline resolutions as done by Stewenius. 

8 Conclusion 

In this paper a new method for the problem of the "five points pose problem" 
has been described. The main difference with the previous methods is that the 
rotation and the translation are directly the unknowns of the system. Then we 
have shown that with the available tools of geometric algebra, this kind of system 
can be solved. The major advantage, when compared to other methods, is that 
it works accurately on all cases of plane scenes, or on couples of images with a 
null base. 
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